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SUMMARY 


Three-dimensional (3-D) finite strips are formulated by combining finite 
element shape functions with beam eigenfunctions. Because of the orthogonality 
of the beam functions, three-dimensional problems are reduced to a series 
of two-dimensional problems, often with stiffness matrices of very narrow 
bandwidth. These require considerably less computer memory and computation 
time to solve. Isoparametric and high order finite element shape functions 
are used in the formulation of the 3-D finite strips. Numerical examples 
such as the static and free vibration analyses of simply supported thick 
plates are presented. Results are compared with existing solutions. Good 
agreements are obtained in all cases. Potential applications of the 3-D 
finite strips include the static and dynamic analyses of voided slabs, thick 
box girders and axisymmetric thick-walled shell structures. 


INTRODUCTION 


Although the finite element method is at present the most powerful 
and versatile numerical approach for structural analysis, the computing 
cost can often be very high. This is particularly true in the case of three- 
dimensional structural analyses. In an attempt to reduce the computational 
requirements of the finite element method, researchers have developed the 
finite strip technique (ref.1), a semi-analytical method that couples simple 
polynomial expressions for one or two directions with beam eigenfunction 
series for the other directions. This reduces a two-dimensional problem 
to one dimension, and a three-dimensional problem to two dimensions. 

Furthermore, because of the orthogonal properties of the eigenfunction series, 
the terms of the series may become uncoupled depending on the type of boundary 
conditions, and the stiffness matrices of each term can be formed, assembled 
and solved separately, resulting in a substantial reduction in computing 
costs. The method is suited for the analysis of structures having regular 
geometric plans and simple boundary conditions, and has been successfully 
applied to the static and the dynamic analyses of slabs, folded plate structures 
and box-girder bridges (ref. 1). 

In this study, 3-D finite strips are formulated in two ways, one by 
coupling isoparametric quadrilateral and triangular plane-stress finite 
element shppe functions with beam eigenfunctions, and the second by coupling 
high order quadrilateral plane-stress finite element shape functions with 
beam eigenfunctions. The stiffness, mass and load matrices are derived 
following standard finite element procedures. Three-dimensional elasticity 
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constitutive equations are used in the derivation of the various stiffness 
matrices. Applications of the 3-D finite strips to some prismatic solids 
such as thick plates are described. Numerical integration using Legendre- 
Gauss or Radau-Gauss quadratures was employed in the derivations. 

SYMBOLS 
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lengths of sides of a quadrilateral 

span of 3-D finite strips 

surface area of 3-D finite strips 

matrix relating strains to displacement components 

matrix containing finite strip displacement functions 

material constant matrix 

individual and assembled consistent load vectors 

vector containing distributed body forces 

thickness of plate 

triangular area co-ordinates 

individual and assembled stiffness matrices 

individual and assembled consistent mass matrices 

finite element shape functions 

number of nodes in finite elements 

vector containing concentrated nodal forces 

point in ? - n space or a concentrated point 

load 

vector containing distributed surface forces 
volume of 3-D finite strip 

displacement components in the x, y and z directions 
Cartesian co-ordinates 

distance along lines of equal n (equation 8) 
individual and assembled nodal displacement vectors 
curvilinear co-ordinates 

finite element rotational degree of freedom =(9v/3x^)^ 
3-D finite strip rotational degree of freedom =(3w/9x^)^ 
mass density 

circular natural frequencies 
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finite element skew symmetric rotational degree 
of freedom =(3v/3x - 3u/3y)/2 

3-D finite strip skew symmetric rotational degree 
of freedom =(3w/3x - 3u/3z)/2 

ISOPARAMETRIC 3-D FINITE STRIPS 


A family of 3-D isoparametric quadrilateral and triangular finite strips 
can be developed by using the plane-stress isoparametric finite element 
shape functions reported by Ergatoudis (ref. 2). Considering only simply 
supported situations in which u=w=3v/3y=0 at the ends, a suitable set 
of displacement functions for a 3-D strip of span a (fig. 1) is 
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The 

most 

simple isoparametric quadrilateral is the four node IPLQ 

quardi 

laterial 

(ref. 

, 2) 

which has linearly varying displacements 

(fig. 2a). 

The 

shape functions for this finite element are simply 
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Other more sophisticated isoparametric quadrilateral elements of the 
same family include the eight node IPQQ quadrilateral (ref. 2) whose dis- 
placements vary quadratically (fig. 2b), and the twelve node IPCQ quadrilateral 
(ref. 2) with cubically varying displacements (fig. 2c). 
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For the Isoparametric triangular finite elements, the shape functions 
are most conveniently expressed in terms of the area co-ordinates L^, L2 and 
L^. The first element of the series is the three node IPCST constant strain 

triangle (fig. 3a) whose shape functions are simply the area co-ordinates 
(ref. 2). Thus, 

= L^, N2 = L2, Ng = (^) 

Using a recurrence formula (ref. 2), more refined triangular elements 
such as the six node IPLST linear strain triangle (fig. 3b) and the ten 
node IPQST quadratic strain triangle (fig. 3c) can be formulated. 

Although in theory more refined elements can be derived by introducing 
additional nodes, such elements are often of limited practical use since 
they usually result in stiffness matrices having very large bandwidths. 


HIGH ORDER QUADRILATERAL 3-D FINITE STRIPS 


Because of its linearly varying displacements, the accuracy of the 
IPLQ 3-D finite strip is usually very limited. The strip could be refined 
by adding extra nodes as was done in the last section. However, such a 
procedure is not always desirable, since the bandwidths of the stiffness 
matrices are increased. The alternative, and probably most effective, 
approach is to introduce additional degrees of freedom at the nodes. Two 
high order plane stress quadrilateral finite elements were selected from 
the published literature for this purpose. 


The first element is the QCC3 in-plane quadrilateral with displacements 
u, V and the skew symmetric rotations 


w 



/ 3v 9u 
^9x “ 3y 


(5) 


as degrees of freedom (fig. 4a). This element was derived by Abu-Ghazaleh 
(ref. 3) and subsequently used by Scordelis (ref. 4) to analyze box-girder 
bridges . 


Using the shape function of this finite element, a high order 3-D finite 
strip with u, v, w and the skew symmetric rotations 
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as degrees of freedom is formulated (fig. 4a). Considering only simply 
supported cases, the displacement functions of the finite strip can be written 
as 
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V = Z 
m=1 

CO 

W = Z 

m=1 

where the shape functions N. . , N„. and N_. are the same as those used for 

lx 2i 3x 

the finite element in reference 3. 

The other high order element selected is the plane stress QLC3 element 
developed by Sisodiya et al. (ref . 5). The nodal parameters of the element 

T 

are {6}. = [u, v, 9 ] ., with 0 .=(9v/9x^). where x_ is a distance along lines of 

1 Z 1 Z1 f, X 5 

equal q (fig. 4b); and at a general point P 

^ {a^ (1-n) + (8) 

where a^ and a^ are the lengths of opposite sides of a quadrilateral (fig. 4b). 
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Adopting the shape function of the element for the 3-D finite strip, 
the displacement functions of a simply supported 3-D strip can be expressed 
as 
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where 9 . = 
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and are 


(9w/9x^) . and x^ is as defined before. The functions 
identical to those presented in reference 5. 
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STIFFNESS, MASS AND LOAD MATRICES 


The stiffness and load matrices can be derived through the minimization 
of the total potential energy, a standard finite element procedure that 
leads to the familar expression 
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[K] {6} - [F] = 0 


( 10 ) 


in which [K] is the stiffness matrix, {5} the unknown nodal displacement 
vector and {F} the consistent load vector. A typical submatrix [K. .] of 
the stiffness matrix [K] is 

iK. .] = /[Bj T[d 1 [B] .dV (11) 

1 J 1 J 

where [B] is the so-called strain matrix that relates the strains to the 
displacement components and [D] is the elasticity matrix for the material 
which can be isotropic or orthotropic. A typical submatrix {F.} of {F} 

is 


{Fi> = {P.} + J[C^]'^{q}dA + /[C^l’’'{g}dV (IP) 

where [C^] contains the nodal displacement functions and the force terms 
represent concentrated, surface and body forces. 


Using the displacement functions defined in the previous sections, 
equation (11) would become 




'D1 [p.; dV 
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Because of the orthogonality of the series used, it can be shown that 
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a ri “ 
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i.e., the series terms are uncoupled and off-diagonal submatrices in the 
stiffness matrix are null matrices. 


To obtain the consistent load vector for the 3-D strips, the external 
applied loads are expressed in terms of series similar to those used for 
the displacement functions and substituted into the appropriate integrals 
in equation (12). Details of the derivation of the load terms can be found 
in reference 1 . 


The formula for deriving the consistent mass matrix is quite standard, 
and is 

[M] = /p[C]^[C]dV (15) 

where [M] is the consistent mass matrix and p is the mass density 
of the material. 

As the displacement functions are either defined in terms of the cur- 
vilinear co-ordinates ^ and q or area co-ordinates L^, and L^, it 

is necessary to rewrite the derivatives and integrals of the displacements 
with respect to the local co-ordinate system. This is a fairly straightforward 
matter involving the determination of the Jacobian matrix (ref. 2). 
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Onoe the individual finite strip stiffness, mass and load matrices 
are formulated, they can be assembled in the usual manner to form the static 
problem of 

[K] {S} - {F} = 0 (16) 

or the free vibration problem of 

([K] - 0)^ [M] ) {S} = 0 (17) 

where [K] , [M] , {f} and {S} are respectively the assembled stiffness, mass, 
load and displacement matrices, and w is the circular frequency of free 
vibration. 


ILLUSTRATIVE EXAMPLES 


Static and free vibration analyses were carried out for the simply 
supported thick, square plate shown in fig. 5. The central deflections 
due to uniformly distributed load and a central point load are shown in 
tables 1 and 2. The. values shown were obtained with ten series terms. 

Comparing the present results with existing finite element solutions (ref. 

6) and closed form solutions (ref. 7) as well as the classical thin plate 
theory (ref. 8), it can be seen that the agreement is good. By doubling 
the number of series terms in a number of runs, it' was found that the displace- 
ment values remained unchanged. The first three lowest flexural frequencies 
for a simply supported square plate with a thickness vs. span ratio of 0.2 
are tabulated in table 3. Excellent agreement between the present frequencies 
and those obtained from finite element (ref. 9) and closed form (ref. 10) 
solutions can be seen, while the thin plate theory tends to overestimate 
the frequencies. The results in tables 1-3 indicate that the effects of 
thickness-shear deformation and rotary inertia can be accurately predicted 
by the present 3-D finite strip formulation. 


CONCLUDING REMARKS 


Three-dimensional (3-D) simply supported finite strips with quadri- 
lateral and triangular cross sections have been formulated using isopara- 
metric and high order finite element shape functions and beam eigenfunctions. 

In general, the accuracies of both the isoparametric and high order 3-D 
strips can be considered good since reasonably good results can be achieved 
even with a relatively coarse mesh. 

Three-dimensional finite strips with other than simply supported boundaries 
can be derived by employing the appropriate beam eigenfunctions to match 
the boundary conditions. Curved, and circular 3-D strips can be developed 
by using a cylindrical co-ordinate formulation. Continuous structures can 
be analyzed either by coupling the finite element shape functions with 
eigenfunctions of continuous beams (ref. 11); or by using a finite strip 
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flexibility approach (ref. 1). Potential applications of the 3-D finite 
strips include the static and dynamic analyses of voided slabs, thick box 
girders and axi symmetric thick-walled shell structures . These potential 
applications along with the above mentipned extension of the 3-D finite 
strip method are the topics of current investigations by the authors; and 
the results will be reported when they become available. 
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TABLE 1 COMPARISON OF CENTRAL DEFLECTIONS FOR 
SIMPLY SUPPORTED PLATES UNDER UDL 



h/a=0.05 

h/a=0. 1 

h/a=0.2 

h/a=0.25 

IPLQ 

3692.1 

463.4 

64.1 

35.4 

IPQQ 

3694.0 

482.1 

65.5 

36.3 

IPCQ 

3704.2 

487.2 

67.1 

36.4 

IPCST 

3681.1 

462.1 

63.2 

34.2 

IPLST 

3688.5 

465.2 

64.8 

36.3 

IPQST 

3701 . 1 

471.3 

66.1 

36.5 

QCC3 

3683.1 

465.4 

64.3 

35.5 

QLC3 

3686.5 

465.1 

64.5 

35.8 

REF. 6 

3575.2 

461.2 

64.8 

35.9 

REF. 7 

3588.8 

463.2 

65.2 

36.2 

REF. 8 

3549.6 

443.7 

55.5 

28.4 


TABLE 2 COMPARISON OF CENTRAL 

DEFLECTIONS FOR 



SIMPLY 

SUPPORTED PLATES UNDER POINT LOAD 



h/a=0.05 

h/a=0.1 

h/a=0.2 

h/a=0.25 

IPLQ 

105.98 

13.97 

2.24 

1.39 

IPQQ 

106.52 

14.47 

2.31 

1.43 

IPCQ 

106.71 

14.87 

2.36 

1.49 

IPCST 

104.92 

13.88 

2.27 

1.39 

IPLST 

106.16 

14.32 

2.29 

1.42 

IPQST 

106.89 

14.77 

2.32 

1.48 

QCC3 

106.11 

14.12 

2.28 

1.41 

QLC3 

106.31 

14.16 

2.30 

1.41 

REF. 6 

106.49 

14.77 

2.46 

1.47 

REF. 8 

101.34 

12.67 

1.58 

0.81 


TABLE 3 COMPARISON OF CIRCULAR 

FREQUENCIES OF A 



SIMPLY 

SUPPORTED THICK 

: PLATE (h/a=0.2) 




MODES OF VIBRATION 



(NUMBEF 

1 OF HALF-WAVES 

IN X AND y DIRECTIONS) 



1,1 

1,2 

2,2 

IPLQ 


1.0621 

2.3324 

3.4000 

IPQQ 


1.0611 

2.3251 

3.2781 

IPLST 


1.0625 

2.3351 

3.3386 

IPQST 


1.0608 

2.3284 

3.2762 

QCC3 


1.0631 

2.3342 

3.4017 

QLC3 


1.0635 

2.3289 

3.3947 

THIN PLATE 

THEORY 

1.1947 

2.9867 

4.7788 

REF. 10 


1.0607 

2.3291 

3.3765 

REF. 9 


1.0565 

2.3235 

3.2758 
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Figure 5.- Mesh sizes us 
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DYWAMIC ANALYSIS- 0=1, E=l,i)=0-3, <0=1 
h 


ESH FDR QUADRILATERAL 3-D STRIPS 



for simply supported thick plate analyses. 



